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Quantum-enhanced measurements use highly non-classical quantum states in 
order to enhance the sensitivity of the measurement of classical quantities, like 
the length of an optical cavity The major goal is to beat the standard quan- 
tum limit (SQL), i.e. a sensitivity of order where N is the number of 
quantum resources (e.g. the number of photons or atoms used), and to achieve 
a scaling known as the Heisenberg limit. Doing so would have tremendous 
impact in many areas , but so far very few experiments have demon- 
strated a slight improvement over the SQL The required quantum 
states are generally difficult to produce, and very prone to decoherence. Here 
we show that decoherence itself may be used as an extremely sensitive probe of 
system properties. This should allow for a new measurement principle with the 
potential to achieve the Heisenberg limit without the need to produce highly 
entangled states. 

Decoherence arises when a quantum system interacts with an environment with many 
uncontrolled degrees of freedom, such as the modes of the electromagnetic field, phonons in 
a solid, or simply a measurement instrument . Decoherence destroys quantum mechanical 
interference, and plays an important role in the transition from quantum to classical mechan- 
ics js]. It becomes extremely fast if the "distance" between the components of a "Schrodinger 
cat" -type superposition of quantum states reaches mesoscopic or even macroscopic propor- 
tions. Universal power laws rule the scaling of the decoherence rates in this regime 
and lead to time scales so small that in fact the founding fathers of quantum mechanics 
postulated an instantaneous collapse of the wave-function during measurement. Only re- 
cently could the CO 



cat" -states 



11 



lapse be time-resolved in experiments with relatively small "Schrodinger 
12l |. However, different superpositions may decohere with very different 
rates. In particular, if the coupling of the quantum system to the environment enjoys a 
certain symmetry, entire decoherence-free subspaces (DFS) may exist, in which superposi- 
tions of states retain their coherence, regardless of the "distance" between the superposed 
states. In essence, the symmetry prevents the environment to distinguish the states, such 
that no information leaks out of the system and the quantum superpositions remain intact. 
DFS have found widespread use in quan tum information theory after their formulation for 



Markovian master equations Q, 14, 1^, experimental demonstration 
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it was realized that quantum computation might be performed inside a DFS 19|] . Given the 



and once 
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reliance of the DFS on a symmetry in the couphng to the environment, it is clear that for a 
large "Schrodinger cat" -type superposition prepared in a DFS, the decoherence rate should 
be extremely sensitive to any changes that modify the symmetry of the coupling. This is the 
basic idea underlying the new measurement principle which we call "Decoherence- Enhanced 
Measurements" (DEM). Two fortunate circumstances make us believe that this idea may 
be turned into something of practical relevance. First, while it may seem that DEMs would 
again require the extremely difficult initial preparation of a highly entangled macroscopic 
state, surprisingly the Heisenberg limit can be reached with a much simpler to prepare prod- 
uct state of pairs of atoms. Second, the required initial "symmetry" of the coupling to the 
environment means nothing more but a degenerate eigenvalue of the Lindblad operators in 



the master eq 



environment 



nation, or, more generally, of the coupling Hamiltonian of the system to the 



14| . Actual symmetries in the system (e.g. spatial symmetries) may lead to 



such degeneracy, but are by no means necessary 201]. The scheme is therefore much more 
general than it may appear at first sight. 

In order to illustrate the concept, we consider two-level atoms or ions {N even, ground 
and excited states |0)j, for atom i, i = 1, . . . ,N) localized in a cavity with one semi- 
refiecting mirror, and resonantly coupled with coupling constants Qi to a single e.m. mode 
of the cavity of frequency u. The reduced density matrix p of the atoms evolves according 
to the master equation in the interaction picture 

|p(t) = /:.[p(t)] + /:,[p(t)], (1) 

where 

r ^ 



^sipit)] = ] + ]) (2) 

1=1 

describes individual spontaneous emission with rate F, while 

^Mt)] = 7 {[J~p{t), J+] + [J^, p{t)J+]) (3) 

models collective decoherence, J_ = ^+ — '^+^ — ('^-'' = 

and Qi = Qi/g where g is the average coupling strength over all atoms coupled to the cavity 

field. The rate 7 = g^/ n (with k, the single photon cavity decay rate) is independent of A^. 



Equation ([T]) is a well-known and experimentally verified 2l|, |22[ master equation w 



lich for 



gi = and in the bad cavity limit F <^ gy N <^ k describes superradiance 



23 



24 



25 



26|. 



Due to the spatial envelope of the e.m. mode in resonance with the atoms, the gi depend 



on the position Xi of the atoms along the cavity axis and on the length L of the cavity (the 
waist of the mode is taken to be much larger than the size of the atomic ensemble), 



^ sm{k^Xi) e ■ d , (4) 



eoV 

where k.^ = tttIx/L, eo denotes the dielectric constant of vacuum, V = LA the mode volume 
(with an effective cross-section A), e the polarization vector of the mode, and d the vector of 
electric dipole transition matrix elements between the states |0)j and taken identical for 
all atoms. Decoherence in this system has been extensively studied, see 20| for a review. The 
initial state po = |^o)(V'o| belongs to the DFS with respect to collective emission, >CJpnL= 0, 
if and only if J.\ipo) = Q. If = 1 for i = 1, . . . , A^, this DFS is well known 20, 27|. It 



contains (^^2) ~ 2^ /VN DP states, including a 2^/^ dimensional subspace <^^J:i{\t-)i, \s)i} 
in which the pair formed by the atoms / and / + N/ 2 can be in a superposition of the triplet 
ground state = \0)i\0)i+n/2 and the singlet \s)i = ^(|0)j|l)/+Ar/2 - |l)/|0)/+Ar/2). For 
arbitrary (ji, \s)i should be replaced by ,^ {gi\0)i\i)i+N/2 - gi+N/2\'^)i\0)i+N/2)- 

^^/\9l\ +\9l+N/2\ 

Consider now the situation where the atoms can be grouped into two sets with N/2 
atoms each and coupling constants Gi in the first set (z G 6*1 = {1 . . . , A^/2}), and G2 in the 
second set (z G 6*2 = {N/2 + 1, . . . , N}). One way of obtaining two coupling constants may 
be to trap the atoms in two two-dimensional lattices perpendicular to the cavity axis (see 
FigH]). Suppose that after preparing the atoms in a DFS state corresponding to the initial 
couplings Cf^ = Cf^/g (/ = 1, 2), the length L of the cavity changes slightly. The coupling 
constants will evolve, Cf' Gj = Gj/ g, and so will the DFS. It is this collective change 
of the coupling constants which can be revealed very sensitively through the decoherence 
it induces as the original state becomes exposed to decoherence. The induced decoherence 
therefore provides for a very precise measurement of the change of the length of the cavity, 
as we shall show now. 

In order to simplify notation we will assume in the following G^'^ = G^2^ (i.e. the atoms 
are located for instance symmetrically with respect to an antinode of the cavity mode, or at 
a distance given by an integer multiple of the wavelength of the mode), but we emphasize 
that everything goes through for different initial couplings, unless otherwise mentioned. 
Assume that an initial pure product state of pairs of atoms is prepared in the initial DFS, 
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FIG. 1: N atoms (or ions) are trapped at fixed positions in two 2D optical lattices perpendicular to 
the cavity axis. A dipole transition of the atoms is in resonance with a single, leaky cavity mode. 
The atoms are initially prepared in a DFS state relative to a given cavity length. When the cavity 
length changes slightly, the DFS evolves, and the initial state is exposed to collective decoherence, 
detectable by photons leaking out through the semi-reflecting mirror at a rate proportional to N'^. 

Po = \i^o){i^o\, where 

N/2 

l^o) = (g) (5) 
1=1 

with \(pi)i = ai\t^)i + bi\s)i. The decoherence mechanism (j3]) is directly linked to photon 
loss from the cavity. The induced decoherence can be measured through the number of 
photons nph which escape through the cavity mirror during a small time interval At. In the 
superradiant regime considered here (F ^ gVN <^ k), any photon created leaves the cavity 
immediately, such that the quantum expectation value of riph is given by (^ph) = —{Jz{t = 
0))cAt = —tr[JzCc[po])At, where the collective pseudo-spin component = | Xlili '^z*'' 
measures total population inversion of the atoms. As long as \iPq) resides in the DFS, we 
have (riph) = 0. If the coupling constants undergo slight changes and get replaced by general 
values gi for atom i, a straightforward calculation (see Methods) shows that 

N/2 

{jz{0))c = -'ji^^lgi - gi+N/2\^\bi\^ 
1=1 

N/2 

+ i^i - gi+N/2) {dj - gj+N/2) b*bjaiafj , (6) 

which is in general of order A^^. In particular, if the Gj undergo collective changes Cf^ —>■ Gj 
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with Gi 7^ G2 and if all pairs of atoms were prepared in the same initial state ai = a, bi = b 
for / = 1, 2, . . . , N/2, we have 

(^.(0))c = -i\5G\'f{N,b) (7) 
f{N,b) = 4|6r[| + |(|-l)(l-|6n] , (8) 

where 5G = {Gi — G'2)/2. The term quadratic in N is maximized for |6| = 1/V2, i.e. an 
equal weight superposition of the two DF basis states and \s)i for each pair of atoms, 
and gives for ^ 1 a signal (?T-ph) — 7At|5G'pA^^/4. As long as the two lattices are not 
situated at anti- nodes of the mode, the relation between 6G and 6L/L is linear to lowest 
order. If we choose X2 — xi = mX with Ux — 1 > m & N we have 

oG = m7rcot( — - — j-p, (9) 

where we see that 5G and 5L/ L are related by a factor independent of A^. Note that the mea- 
surement of (riph) allows the measurement of and not just a detection of a change of L : 

/ \ 1/2 

5L/L= ( (?T,ph)/(7At(m7rcot(na;7rxi/L))^/(A^, 6)) ) . The ultimate sensitivity achievable 



depends not only on the scaling of the signal (nph) with A^, but also of the noise, quantified 
through the standard deviation aiji^T^) = {{n'^^) — (?^ph)^)^^^. The most fundamental noise 
associated with the measurement of riph is its fluctuation due to the quantum mechanical 
nature of the prepared state. Our approach of calculating initial time-derivatives of observ- 
ables by tracing them over with the Lindbladian implies (?^ph) ^ 1, as J_ can change the 
number of excitations by at most 1, and this condition sets an upper bound on At. In this 
regime, (^iph) ^ (?^ph), and, therefore, a(?7,ph) ~ \/ (nph). It follows that the signal-to-noise 
ratio is given by (7T.ph)/o"(raph) — \6G\^/^yAt^^/J{NJ)}. The ultimate sensitivity achievable 
can be estimated from a fixed (riph)/o"('^ph) of order 1, independent of A^, which leads to a 
minimal \SG\ = 0{l/{y/'~fAtN)). We have thus shown that a precision measurement based 
on the purely dissipative dynamics ([3]) and an initial product state can achieve the Heisen- 
berg limit. This is in contrast to unitary dynamics of A^ independent quantum resources. 



where the SQL cannot be surpassed when using an initial product state [28|. 

In a real experiment there may be additional fundamental noise sources. One obvious 
concern is spontaneous emission. It is easily verified that £s[po] leads to a contribution 
{jz{0))s = ~'^'l2f=i l^iP (^^(0)) "which scales as 0{N) and leads to a background signal 
against which, one might think, the collective decoherence signal {Jz{0))c has to be com- 



pared. However, note that spontaneous emission sends photons into the entire open space 
but not into the cavity, whereas the collective emission escapes exclusively through the leaky 
cavity mirror. Therefore, the two contributions can be well separated by observing only the 
photons escaping through the cavity mirror. 

Another obvious concern are fluctuations of the coupling constants. In order to measure 
(riph), the experiment has to be repeated u times with u ^ 1. However, u is independent 
of and only given by the desired signal/noise, or, equivalently, by (?^ph) itself. Increasing 
u does therefore not influence the scaling with A^. It increases the sensitivity by a factor 
V^, but reduces the bandwidth by 1/z/ in the standard way. While we assume that the time 
scale of the mirror motion is sufficiently long compared to the time needed for averaging, 
the exact coupling constants might fluctuate about their slowly evolving mean values during 
the averaging, e.g. due to fluctuating traps caused by vibrations in the set up. But even for 
perfectly stable traps, the micro-motion of the atoms in their respective trapping potentials, 
thermal motion, or even quantum fluctuations in the traps will lead to fluctuating gi. The 
cost in sensitivity of these fluctuations depends on their correlations. To see this, let us 
consider fluctuations 6gi of the Qi about their mean values Gj, Qi = Gi + 6g.i for i = 
1, . . . , N/2, Qi = G2 + Scji for i = N/2 + 1, . . . , A^. We introduce the correlation matrix 
Cij = 6g*6gj, where the over-line denotes an average over the ensemble describing the 
fluctuations. Equation ([6]) then leads to a background abg in the photon counting rate 
and to fluctuations 6af on top of ac = —{Jz{0))c, i.e. a = ac + «bg + <^a/- The average 
background abg, given by 

N/2 
1=1 

can be determined independently at 6G = 0, and subtracted from the signal; it does therefore 
not influence the sensitivity of the measurement. The remaining noise 6a f fluctuates about 
zero, 

N/2 N/2 

Saf = 27{ J2i^G*A~g, + 6GA~g*)\k\' + J2i^G*A~g, + 6GAg*)b*bja,a*} , (10) 

i=l *,j=i 

where Acji = 5gi — 5giJ^N/2- Assuming real coupling constants, we find the standard deviation 
(5^)1/2 = A^\5G\K, where K = (Et^=i(C,, + Q+iv/2j+^/2 - Q,,+^/2 - C,^N/2,)S^S,f'\ 
and Si = + biO'iJ2j\jj!=ibj'^j general of order 0{N). This leads to fiuctuations in 
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the measured photon number with a standard deviation (Sn^^Y^'^, and to a signal-to- noise 
ratio (nph)/(5^)^/2 = |5G|/(iV, 6)/(4K). 
Several interesting cases can be considered: 

1. Fully uncorrelated fluctuations, Cij = CiSij, where Sij stands for the Kronecker-delta: 

1 /2 

Here we get K = (j2f^i{Ci + Ci+N/2)Sf^ , which is in general of order N'^^"^, and 
leads back to the SQL. 

2. Pairwise identical fluctuations between the two sets: Ca = C-,n_- = C--,n_ = 
C^_^Nj^N for i,j = 1, . . . , N/2. This can be the consequence of fully correlated fiuc- 
tucations, Cij = C V Alternatively, such a situation arises for example for atoms 
initially arranged symmetrically with respect to an anti-node such that = 0^2^ 
if the two atoms (or ions) in each pair / (/ = 1, . . . , N/2) are locked into a common 
oscillation. This should be the case for two trapped ions repelling each other through 
a strong Coulomb interaction, and cooled below the temperature corresponding to the 
frequency of the breathing mode. Equation (fTOl) then gives 6af = 0. Note, however, 
that for initial cf^ ^ Cf the more general DFS leads to a more complicated condition 
for the correlations, Cij\Gf\'^+C^^N.j^N\Gf\'^-Cij^NGf*Gf^-Ci^ = 0, 
which might be harder to achieve. 

3. Correlated fluctuations within a set, but uncorrelated between the two sets, Gij = G 
for i,j G 5*1 or i,j G 5*2, but Gij = for z G 5*1 and j G 5*2 or vice versa. In this case 
both sums in (ITOl) survive and lead to a noise of order (9(A^^), the worst case scenario. 
However, this comes to no surprise, as such correlations are indistinguishable from the 
signal: all the atoms in a given set move in a correlated fashion, but independently 
from the atoms of the other set. This leads to a collective difference in the couplings, 
just as if the length of the cavity was changed. 

Case (2) above is clearly the most favorable situation. If there are no other background 
signals depending on A^, we keep the scaling of 6G characteristic of the Heisenberg 
limit. In order to favor case (2) over cases (1),(3), it appears to be advantageous to work 
with ions and to try to bring the ions in a pair as closely together as possible, thus strongly 
correlating their fluctuations, while separating the ions in the same set as far as possible. 

To summarize, we have shown for a particular example how the very sensitive dependence 
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of collective decoherence on system parameters can be exploited to reach the Heisenberg limit 
in precision measurements while using an initial product state — something which is known 
to be impossible with unitary dynamics 28|]. It should be clear that the principle of DEM 
is far more general than the example exposed here. Decoherence is itself a process in which 
interference effects play an important role. This is exemplified by the very existence of DFS, 
and can lead to exquisite sensitivity. One might therefore as well try to exploit these effects 
instead of trying to suppress decoherence at all costs. 
Methods 

Derivation of Eq. (Q: The commutation relation [J±,Jz] = 1^J±, valid for any choice of 
couphngs cji, allows to rewrite {jz{0))c = ~2j{iPq\J+J_\iPq). A short calculation yields 
J-l^po) = -{'^/V^)Yli'=ii9i - 9i+f)bi\t-)i<S)'i^-li and leads immediately to Eq. 
Preparation of initial state: In order to prepare the product state ([5]) it is helpful to use 
three-level atoms with a lambda structure. Let |0) and the additional state |2) be hyperfine 
(HP) states, and assume that their energies are split in a sufficiently strong magnetic field, 
such that only the transition |0) ^ |1) resonates with the cavity mode. We assume further 
that the second optical lattice can be moved along the cavity axis, such that controlled 
pairwise collisions of corresponding atoms in the two lattices can be induced. Entangled 
pairs of atoms in their HP split ground states can thus be created (for atoms in the same 



lattice this has been demonstrated experimentally, see [29| for a review). After the creation 
of an entangled HP state \iPq), that differs from ([5]) by the replacement of states |1) by states 
1 2), the second lattice is moved back to its original position. Now one can selectively excite 
the |2) states by a laser pulse in resonance with the |2) ^ |1) transition, that replaces the 
singlets in the (very long lived) HP states by the desired singlets of the |0) and |1) states 
and thus produce the product state ([5]). However, as such, the method is not of much 
practical use yet, as it will be virtually impossible to park the second lattice at the exact 
position corresponding to coupling constants which render the state ([5]) decoherence free. 
The extreme sensitivity of the collective decoherence with respect to changes of the coupling 
constants plays against us here, and will lead to a superradiant flash of light from the cavity 
after the excitation \iPq) I'lpo), if the exact position corresponding to \ipo) G DPS is not 
achieved. But it is possible to position the second lattice at the required position with a 
precision of 0{1/N) for the case (2) considered above, using a feed-back mechanism and a 
part of the quantum ressources. With the atoms in the state \iPq), do the following repeatedly 
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in order to find the optimal position: Excite a part of tlie entangled HF pairs containing 
0{N/ In N) atoms with the laser, measure (riph), and use the measurement results to bracket 
the minimum of (riph) as function of the lattice position. The minimum of (^ph) indicates 
that the position corresponding to the DFS is achieved. Using golden section search, the 
minimum can be bracketed to precision in 0(lnN) moves, as at each step the sensitivity 
of the measurement of the position of the lattice is of order 0{lnN/N) ~ 0{1/N). Once 
the minimum is found, excite the remaining unused pairs (there should be still a number 
of pairs of 0{N)) to the desired state \ipo). That state is now decoherence-free, and the 
system ready to detect small changes of the position of one of the mirrors. Note that for this 
method it is not necessary to know which exact state is produced in the controlled collisions 
and subsequent laser excitation. 

Imperfections in preparation of {ipo): Suppose that instead of the state ([5]) a state 

N/2 

li'o) = 6d with = a\t_)i + b\s)i + c\to)i + d\U)i (11) 



1=1 

was prepared (we consider the same state for all pairs for simplicity, but this is not essential). 
Repeating the calculation that leads to Eq. ([H]) and assuming real G^, we now find 

(riph) = ^M[^[{{c-h)G, + {c + h)G2Y + d\Gl + Gl)) 

+ y - ((G2 - Gi)h{a - d) + (G2 + Gi)c(a + d))'] . (12) 

The derivative of (wph) with respect to G2 is of order 0{N^), and thus still allows to find 
the minimum of a as function of the position of the second lattice with a precision of order 
0{1/N) in case 2. At the minimum a component outside the DFS persists, such that 
photons will leak out of the cavity, but the average background is only of order 0{N), and 
can be measured separately and subtracted from the singal. Changes 56*2 of G2 away from 
the position of the minimum still lead to a signal that scales, for large A^, quadratically 
with A^, (riph) = ^ (a(c + b) + d{c — h))^ N'^5G\, and the previous analysis leading to the 
Heisenberg limit still applies. 

Alternatively, one can get rid of the additional background by letting the system relax 
before measuring changes of L. Indeed, any state with a component outside the DFS will 
relax to a DFS state (and thus a dark state) or mixtures of DFS states within a time of 
order I/7 or less. Components with large total pseudo-angular momentum J relax in fact 
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in much shorter time of order 1/(J7). The DFS states reached through relaxation starting 
from still allow a scaling of {n^h) close to A^^. We have shown this by simulating the 
relaxation process with the help of the stochastic Schrodinger equation (SSE) corresponding 
to Eq. For real values of the coupling constants, the SSE reads 

di){t) = Di{tp{t)) dt + D2{tlj{t)) dW{t) (13) 
Di(^) = 7 (2(J_)^J_ - J+J- - (J_)y ^ (14) 
D^i^) = v^(J_-(J_)^)^/;, (15) 

where dW{t) is a Wiener process with average zero and variance dt, and (J-)^ = {ip\J^\ip) 
jsol . Using an Euler scheme with a time step of 0.01 / 7, we followed the convergence of tlj{t) to 
DFS states for states with (a, b, c, d) = (cos 5, 1, 0, sin S) /V2, until the norm of the difference 
\ip{t + dt)) — {ipit)) dropped below 10~^^. In these final dark states, randomly distributed 
over the DFS, we calculated a, and averaged over a large number rir of realizations of 
the stochastic process (n,. = 10^ 10^ 10^ 2.5 ■ 10^, 10^, 1.25 ■ 10^, 250, 200, and 250 for 

= 2, 4, 6, 8, 10, 12, 14, 16, and 18). Figure [2] shows the scaling of a as function of for 
different values of 5 for < 5 < 7r/2 up to = 18. Within this numerically accessible range 
of A^, a follows a power law a oc A^^ with an exponent p that decays only gradually with S 
for 5 < 7r/4. Moreover, that decay might be a finite size effect: Note that, surprisingly, a{6) 
appears to be close to symmetric with respect to 6 = ir/A. This is corroborated by exact 
analytical calculations based on the diagonalization of Cc, which lead to a = 1/2 for A^ = 2, 
a = (55 - 12sin(25) - cos(45))/36 for A^ = 4, and a = (303 - 110sin(25) - 3 cos(45))/100 
for A^ = 6 (in units 716*1 — 6*2 P). The plot shows that all numerical data can be very well 
fitted hj a = A + B sin(2(5) + C cos{46). From Eq. ([8]) we know that A + C has to scale as 
A^^ for sufficiently large A^. Both B and C are negative for all A^ for which we have data, 
and C appears to be negligible. Fig. [2] shows that —B increases even more rapidly than A^^ 
(a fit in the range A^ = 8, . . . , 18 gives a power law N^'^). But B has to cross over to a 
power law A^-*^ with p < 2, unless other Fourier components start contributing significantly. 
Otherwise, a would become negative for 6 > 0. This indicates that for large A^ the scaling 
of a is in fact A^^ for all 6. 

In summary, our method still works, even if the product state ([5]) is not prepared perfectly. 
One has the choice to start measurement immediately after state preparation, which gives 
an additional background of order A^, or to wait a time of the order of a few I/7 after 
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FIG. 2: (a) Scaling of collective photon emission rate a = —{jz{0))c (in units 'y\Gi — (52p) in the 
mixture of DFS states reached by relaxation from state |^o)) Eq- (|11P - Data for 6 = 5o (5 = 7r/2— Jq) 
denoted by squares and full lines (diamonds and dashed lines), respectively; = 0, vr/9, 27r/9, 37r/9 
and 47r/9 in black, red, green, blue, and brown. Exact analytical results for 6 = 0, Eq. ([8|) shown 
with black crosses. Full and dashed lines are guides to the eye only, (b) Dependence of a on 6. 
N = 2, 4,6,8,10,12,14,16,18 in black circles, red squares, green diamonds, blue triangles up, brown 
triangles left, grey triangles down, violet triangles right, cyan pluses, orange Xs, respectively. The 
full lines for A^=2,4,6 are exact analytical results. The dashed lines for = 8, . . . , 18 are fits to 
A + B sm{25) + C cos(4(^). The inset shows the scaling of the coefficients A (black circles), —B (red 
squares), and —C (green diamonds) as function of N. 



preparation of the initial state, until no more photons leave the cavity through the mirror, 
with no additional background. In both cases the scaling of the signal-to-noise ratio is 
still ~ y/''yAt\6L/L\N, and allows to reach the Heisenberg limit in the measurement of a 
subsequent small change of L. 

Another class of states in the DFS that allows quadratic scaling of a with A^, are 
Schrodinger cat states of macroscopic pseudo-angular momentum ji = j2 = i ^ N in 
the two sublattices (i.e. states — j)), 



-j\Ji+Ji-\{i,i)j, -J) = \G,\ 



,(2j + l)(2j + 2) 



1) 



j(j+2) 



(16) 



4(j + 1)2-1 

with Ji_ = Gi Xlil^i i'^ agreement with the initial intuitive reasoning. If the total system 
is initially in a singlet state (j = 0) and the angular momentum of each of the two sets of 
N/2 atoms has its maximal value i = N/A, we have 

'J At. 



■\5G\' N{N + 4:) . 



(17) 
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While these states are protected by the DFS and thus do not suffer the fate of rapid decoher- 
ence of the highly entangled states proposed for QEM, it appears to be still more challenging 
to produce them compared to the product states Our numerical simulations also show 
that states chosen randomly inside the DFS lead on the average only to (?T-ph) oc A^. There- 
fore, it is rather remarkable that the product states and the above Schrodinger cat states 
share the property of scaling of (nph) oc A^^. 
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